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I.  SUMMARY  OF  WORK  AND  DESCRIPTION  OF  WORK  ACCOMPLISHED 


This  instrumentation  grant  supported  the  purchase  of  two  multiprocessor,  high¬ 
speed  UNIX  workstations  (an  HP  J5600  and  an  IBM  RS6000/44P).  These  computer 
system  were  (and  continue  to  be)  employed  in  the  theoretical  study  of  the  interactions  of 
boron  and  aluminum  atoms  with  molecular  hydrogen.  An  understanding  of  these 
interactions  is  of  great  importance  for  the  ultimate  development  of  novel  propellants 
based  on  the  doping  of  light  atoms  into  cryogenic  hydrogen.  Quantum  chemical 
calculations  of  the  relevant  potential  energy  surfaces  were  carried  out.  With  these 
potential  energy  surfaces,  the  energetics,  structure,  and  nuclear  dynamics  of  the  weakly 
bound  clusters  was  investigated,  in  particular  with  fully-quantum  determination  of  the 
bound  state  energies  and  wavefunctions  of  the  binary  complex  and  with  path-integral, 
quantum  molecular  dynamics  simulations  of  more  extended  systems. 

Our  AFOSR  funded  research  project  involves  close  collaboration  with  the 
experimental  group  of  Dagdigian  at  The  Johns  Hopkins  University.  They  have  used 
optical  spectroscopy  to  characterize  experimentally  the  non-bonding  interactions  of  the 
aluminum  atom,  in  both  its  ground  and  excited  electronic  states,  with  molecular 
hydrogen. Recently,  several  excited  electronic  transitions  in  the  binary  Al— H2  complex 
were  observed  and  characterized.^  The  rotational  structure  of  several  bands  were 
resolved  and  assigned,  and  the  average  AI-H2  separation  in  the  ground  vibronic  level  was 
determined  in  order  to  provide  direct  information  on  the  length  scale  of  the  AI-H2 
interaction.  Using  ground-state  combination  differences  and  employing  optical-optical 
double  resonance  spectroscopy,  Dagdigian  and  Tan  were  able  to  determine  the  rotational 
constant  B  and  parity  doubling  parameter  p  for  the  ground  bend-stretch  levels  of  Al(3/i)- 
0H2  and  Al(3p)-/?D2.  These  are  the  most  strongly  bound  nuclear  spin  modifications  of 
the  complexes."^ 

To  understand  and  interpret  these  experiments,  we  used  fully  quantum 
calculations  of  the  ground  states  of  these  binary  complexes.  Two  sets  of  potential  energy 
surfaces  were  employed:  those  computed  earlier  Williams  and  Alexander,^  for  which  the 
H-H  bond  distance  r  was  fixed  at  the  diatomic  equilibrium  separation,  and  those 
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determined  by  new  calculations  in  which  the  dependence  on  r  was  included.  These  ab 
initio  calculations  were  carried  out  with  the  MOLPRO  2000  program  suite, ^  installed  on 
our  new  Unix  workstations.  The  determination  of  the  bound  states  was  done  with  our 
Hibridon  4.1  program  suite,^  developed  partially  with  AFOSR  support. 

The  agreement  of  the  experimental  values  of  the  parameters  and  those  calculated 
on  the  basis  of  PESs  with  the  r  dependence  included  is  excellent,  considerably  improved 
over  the  values  determined  from  the  fixed-r  PESs.^  In  addition,  the  dissociation  energy 
of  the  complex  computed  with  the  fixed-r  PESs  is  -15%  smaller  than  the  value  computed 
with  the  full  PESs.  In  the  AI-H2  complex,  which  has  a  T-shaped  equilibrium  geometry, 
the  increased  attraction  in  the  expansion  phase  of  the  H2  vibration  more  than  compensates 
for  the  reduced  attraction  in  the  compression  phase.  One  conclusion  of  this  comparison  is 
that,  for  weakly  bound  complexes  of  atoms  with  molecular  hydrogen,  it  is  necessary  to 
take  into  account  the  dependence  of  the  PESs  on  the  H2  bond  distance. 

In  a  collaborative  project  with  Soonmin  Jang  and  Greg  Voth  at  the  University  of 
Utah,  we  incorporated  our  B— H2  PESs  into  their  path-integral  Monte-Carlo  simulations  of 
the  diffusion  and  recombination  of  B  impurities  in  solid pYLi-  This  PES,  which 
incorporates  the  dependence  of  the  interaction  potential  on  the  orientation  of  the  2p 
electron,  was  based  on  the  ab  initio  calculations  of  Williams  and  Alexander,'^  is  a 
considerable  improvement  over  the  simplified  potential,  obtained  by  averaging  over  all 
directions  of  the  B  2p  orbital,  which  was  used  in  the  earlier  work  of  Jang,  Jang,  and 
Voth.7 

We  found  that  the  inherent  anisotropy  of  the  hexagonal  close-packed  lattice  is 
significant  enough  to  warrant  the  use  of  electronic-orientation-dependent  potentials 
whenever  possible;  the  average  B— H2  attractive  potential  energy  is  lowered  by  nearly 
20%  when  this  is  included.  When  additional  anisotropy  is  introduced  by  the  presence  of 
a  site  vacancy,  the  importance  of  the  orientation-dependent  potential  is  further  evidenced 
by  the  behavior  of  the  B-/)H2  pair  correlation  functions.  The  orientation-dependent 
potential  introduces  additional  flexibility,  which  allows  the  electronic  and  spatial 
coordinates  of  the  B  to  adapt  to  the  presence  of  an  adjacent  vacancy. 
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Similar  simulations  for  embedded  A1  impurities  in  solid  H2  are  now  underway , 
based  on  our  latest  AI-H2  PESs.^  As  described  earlier  in  this  report,  these  PESs  include 
an  averaging  over  the  vibrational  motion  of  the  H2  moiety,  and  are  considerably  more 
accurate  than  the  earlier  PESs  of  Williams  and  Alexander,^  which  were  determined  for  a 
rigid  H2  distance. 
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